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ABSTRACT 

Investigation of propagation, conversion, and scattering of MHD waves in the 
Sun is very important for understanding the mechanisms of observed oscillations 
and waves in sunspots and active regions. We have developed 3D linear MHD 
numerical model to investigate influence of the magnetic field on excitation and 
properties of the MHD waves. The results show that the magnetic field can sub- 
stantially change the properties of the surface gravity waves (/-mode), but their 
influence on the acoustic-type waves (p-modes) is rather moderate. Comparison 
our simulations with the time-distance helioseismology results from SOHO/MDI 
shows that the travel time variations caused by the inclined magnetic field do 
not exceed 25% of the observed amplitude even for strong fields of 1400-1900 
G. This can be an indication that other effects (e.g. background flows and non- 
uniform distribution of magnetic field) can contribute to the observed travel time 
variations. The travel time variations caused by the wave interaction with mag- 
netic field are in phase with the observations for strong fields of 1400-1900 G 
if Doppler velocities are taken at the height of 300 km above the photosphere 
where plasma parameter f3 <^1. The simulations show that the travel times only 
weakly depend on the height of velocity observation. For the photospheric level 
the travel times are systematically smaller on approximately 0.12 min then for 
the hight of 300 km above the photosphere for all studied ranges of the magnetic 
field strength and inclination angles. The numerical MHD wave modeling and 
new data from the HMI instrument of the Solar Dynamics Observatory will sub- 
stantially advance our knowledge of the wave interaction with strong magnetic 
fields on the Sun and improve the local helioseismology diagnostics. 
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Introduction 



Local helioseismology has provided important results about the structures and dynamics 



regions (e.g. 


Duvall et al. 


1996; 


Kosovichev 


1996; 


Kosovichev et al. 


2000 


Haber et al. 


2OO0I: 


Komm et al. 


2008). The helioseismic inferences help 



the complicated processes of the origin solar magnetic structures, formation and evolution 
of sunspots and active regions. These studies are based on measurements and inversions 
of variations of acoustic travel times and oscillation frequencies in the areas occupied by 
magnetic field and around them. 

There are several factors that may cause the observed variations of the oscillation prop- 
erties, and it is very important to pai nstakingly inv estigate their effects for improving the 
reliability of the helioseismic inference flBogdanll2000l ). Such studies are carried out both ob- 
servationally by doing various experiments with the data analysis proced ure, e.g. by masking 



the r egions of strong field, doing "double-skip" experiments etc. (e.g. iZhao fc Kosovichev 
20061 ). and theoretically by simulating wave propagation in various conditions of the solar 
convection zone and calculating how these conditions affect the helioseismic obseryables, 
such as the oscil lation power spectrum and acoustic travel times (iGeorgobiani et al.l 120071 ; 
Zhao et al.l 120071 ). We emphasize that for correct interpretation of helioseismic results the 
theoretical modeling must include calculations of the actual observables, taking into account 
all important aspects of the data measurement pro cedure, such as data filtering and av- 
eraging, and geometrical factors flNigam et al.l 120071 ). Unfortunately, in many theoretical 
studies the actual helioseismic measurement procedure is not modeled, and this may lead to 
incorrect conclusions about the role of various factors in the helioseismic results. 

In general, the main factors causing variations in helioseismic travel times in solar mag- 
netic regions, can be divided in two types: direct and indirect. The direct effects are due 
to the additional magnetic restoring force, which changes the wave speed and may trans- 
form acoustic waves into different types of MHD waves. The indirect effects are due to the 
changes in the convective and thermodynamic properties in magnetic regions. These in- 
clude depth-dependent variations of temperature and density, large-scale flows and changes 
in wave source distribution and strength. Both direct and indirect effects may be present in 
the observed travel-time and frequency variations and cannot be easily disentangled by data 
analyses causing confusions and misinterpretations. Thus, it is very important to investigate 
the various aspects by numerically modeling the individual factors separately. 

In particular, we have investigated the effects of the suppressed excitation of acoustic 
waves in sunspot regions, where s trong magnetic field inhibits conv ective motions, which are 
the primary source of solar waves (jParchevsky &: Kosovichevll2007a] ). The results showed that 



3 



the suppression of acoustic sources may explain most of the observed deficit of acoustic power 
in sunspot regions (jParchevsky fc Kosovichevll2007bl ). and also cause systematic shifts in the 
travel-time measurements. However, these shifts are significantly s maller than the observe d 
variations of the travel times and and also have the opposite sign (jParchevsky et al.ll2008l ). 
and, thus, cannot affect the basic conclusions on th e sunspot sound-speed structure, contrary 
to previous suggestions (e.g. iRajaguru et al.ll2006l ). 



The goal of this paper is to model the excitation and propagation of helioseismic waves 
(both f- and p-modes) in the presence of inclined magnetic field and investigate the impor- 
tance of the inclined field in the time-distance helioseismology measurements by carefully 
modeling the measurement procedu re and comparing with the observational results, ob- 
tained by IZhao fc KosovichevI (|200(: ]_. The issue of the influence of the inclined magnetic 

field was raised by ISchunker et al.l (120051 ) . who found that the phase shift of the signal in 
the penumbra of a sunspot, measured by the acoustic holography technique varies with the 
sunspot position on the disk. They attributed this to the variations of the angle between 
the inclined magnetic field of the penumbra and the line-of-sight. They suggested that the 
variations of the phase shift may affect the inferences of the sou nd-speed distribution below 
sunspots, inferred by time-distance helioseismology. However, I Zhao &: KosovichevI (120061 ) 
repeated the analysis of the same sunspot by the time-distance technique, and found sub- 
stantially smaller variations with the position on the disk, and no significant effect on the 
wave-speed profile. They also found that the variations due to the inclination angle exist 
only for the wave measurements using the Doppler-shift signal, and that the variations are 
absent when the travel times are measured from the simultaneous intensity observations from 
SOHO/MDI. This result indicates that the observed variations with the inclination angle of 
the Doppler-shift measurements are likely to be related to changes of the ratio between the 
vertical and horizontal components of the displacement vector of the solar oscillations in the 
penumbra, and not the wave transformation or other effects, which could affect the modal 
structure of the oscillations. The solar oscillation theory predicts that the ratio between 
the vertical and horizo ntal components mostly depends on the surface boundary conditions 
(e.g. lUnno et al.lll989l ). In the sunspot umbra, the boundary conditions may change due to 



the inclined magnetic field or/and near surface flows, the Evershed effect, which is observed 
directly in the Doppler-shift data and shows a significant center-to-limb variation. 

In this paper, we present the results of numerical modeling of the inclined magnetic 
field on the time-distance helioseismology measurements by isolating this affect in a sim- 
ple magnetic config u ration , and show that only 25% of travel time variations measured by 
Zhao fc KosovichevI (120061 ) can be explained by a direct influence of the inclined magnetic 
field on acoustic waves. In Sec. 2, we present the governing equation and describe the 
numerical method of 3D modeling of helioseismic MHD waves. In Sec. 3, we present the 
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code verification results comparing the numerical results with analytical solutions for sim- 
ple cases. In Sec. 4, we present the simulation results of the wave propagation in regions 
with inclined magnetic field, calculation of the center-to-limb variations of the time-distance 
helioseismology measurements and comparison with the observational results. 



Numerical Model 



2.1. Governing equations and numerical scheme 

Propagation of MHD waves inside the Sun in the presence of magnetic field is described 
by the following system of linearized MHD equations: 



dp' 
dm' 



V ■ m' = 0, 



4vrpo 



[(V X Bo) xB' + {Vx B') X So] = g^p' + 



dp' 2^ 



so 



Ml 
9q 



0, 



where m' = pov' is the momentum perturbation, v', p', p', and B' are the velocity, density, 
pressure, and magnetic field perturbations respectively, S is the wave source function. The 
quantities with subscript 0, such as gravity g^, sound speed c^o, and Briint-Vaisala frequency 
A/q correspond to the background model. The spatial and temporal behavior of the wave 
source is modeled by function f{x, y, z, t): 



f{x,y,z,t) 




.2 n2 



(l-2r^)e-"' iir<R, 



(2) 



where Rsrc is the source radius, r = ^/ {x — Xsrc)"^ + {y — Vsrc)'^ + ~ ^src)'^ is the distance 
from the source center, r is given by equation 



r 



Lj(t - to) 



vr. 



to<t<to + 



47r 



(3) 



where u is the central source frequency, to is the moment of the source initiation. In our 
simulations we used sources of two types: source of vertical force S = (0, 0, /)^ and pressure 
source S = V/. Superposition of su ch randomly distributed sources d escribes very well the 
observed solar oscillation spectrum (jParchevsky fc Kosovichevll2007al ). 



- 5 - 



For numerical solution of equations ([T]) a semi-discrete finite difference scheme of high 
order was used. At the top and bottom boundaries non-refiective boundary conditions based 
on the perfectly matched layer (PML ) technique were set. We used the standard solar model 
S (IChristensen-Dalsgaard et al.lll996l ) with a smoothly joined model of the chromosphere of 
Vernazza et al. ( 19761 ). as the background model. The background model was modified near 



the photosphere to make it convectively st able. Details of numerical r e alizati on of the code 
and the background model can be found in iParchevsky fc KosovichevI (l2007al ). 



2.2. Code verification for different types of MHD waves 

To verify the code and estimate the total error of the method we compare numerical 
results with simple analytical solutions. We assume that all quantities depend on time t and 
one spatial coordinate x only, and consider the following linearized adiabatic ID system of 
the MHD equations in Cartesian coordinates for a uniform background model 

dp du ^ 
dt ^ dx 

du ^1 dp ^ Bz dBz ^ 



dt p dx Anp dz 
dv Br dB„ 



dt Anp dx 
dw B,. dB^ 



0, 



(4) 



dt Airp dx 
9By _ ^^dv ^ 
dt ^ dx 

^ + Bz^-B^^ = 0, 
dt dx dx 

d-t^'^Pdi = '^ 

where Cg is the sound speed in the background model, p, p are the pressure and density, u, 
V, and w are the velocity components, B^, By, B^ are the components of the magnetic field 
respectively. The quantities related to the background model are marked here by overbar. 
By rotating of the coordinate frame around OZ axis we set By = 0. Equations (jlj) can be 
rewritten in matrix notations 

dU dU 

dt dx 
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where U = {p,u,v,w, By, Bz,p)'^ . We seek a solution of equations ([5]) in infinite interval 
— oo < X < oo in the form of plane waves 

U = Uexp[i{kx - ut)] (6) 

Substituting equation (Q into system ([5]) we obtain an eigenvalue problem for amplitude C/ 

AU = ^U, (7) 

where V = u/k is the phase speed of the wave. Thus, the amplitudes of various MHD 
quantities are the components of eigenvectors of matrix A, and eigenvalues of this matrix 
represent phase velocities of the corresponded waves. From this we calculate the amplitudes 
of the entropy, Alfven, slow, and fast MHD waves: 

Us = ip, 0, 0, 0, 
L/^A = (0, 0, ±Va, 0, 

, ^ c^Vssin6'cos6' 

U^s = ip, TVs, 0, T ^ , 

A S 

, ^ ciVp sin 6 cos 9 

U^F = {p, tVf, 0, T^^-rF2 — 772 — ' 

The phase velocities of these waves are 0, =fVa, =fKs) and T-Vp respectively, where 

Va = cos6', 
1 



0, 


0, 




or 




0, 




or 


0, 




sin Q 




VI- 




0, 




sin Q 




VI- 





Vs = ^\jcl + c\-sjct + c\- 2clc\ cos 29, 
1 

V2 



Vf = —\Ic^ + cl + Jd + 4 - 2c2ci cos 2^. 



Here ca = B / y/^np is the Alfven speed, 9 is the angle between the wave vector and the 
background magnetic field {B^ = Bcos9, B^ = Bsin9). Equations ([6]), ([8]), and Qj give us 
an analytical solution to ID system of MHD equations (jl]), which includes all types of MHD 
waves. Initial conditions for these waves are obtained by setting t = in equation ([6]). 

For testing the numerical simulations we use our 3D code with initial conditions de- 
pending only on x variable. The calculations are carried out in the Cartesian geome- 
try in the domain of 15.46 Mmx 15.46 Mmx3.05 Mm with the numbers of grid points: 
Nx = Ny = 104, Nz = 71. All boundary conditions are chosen periodic simulating an 
infinite spatial domain. Wave vector k = 107r/(A'a; + l)Aa: is chosen in a way to match the 
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periodic boundary conditions. In dimensionless variables p = p/ P-, P = p/ pc^si v = v/cs, and 
= B'^/Airpc^ parameters of the background model are 

p = 1, Cs = l, B = 0.5, ca = 0.5, 9 = 7r/4. (10) 

The amplitude of the Alfven wave (in chosen dimensionless variables) traveling in the positive 
direction of the a;-axis is U+A = (0, 0, -1, 0, 1, 0, 0)^. For the slow and fast MHD waves trav- 
eling in the same direction these amplitudes are U+s = (1, 0.33108, 0, 2.6894, 0, —2.5184, 1)^ 
and tZ+i? = (1, 1.0679, 0, —0.13146, 0, 0.39708, 1)"^ respectively. The cyclic frequencies of 
these waves are uja = kVA, ujs = kVs, and up = kVp respectively. The results of our nu- 
merical and analytical solutions for the moment of time t = 20 min. for all three waves are 
shown in Figure [T] Panels a, b, and c represent the results for the Alfven, slow, and fast 
MHD waves respectively. Only variables v and By are nonzero in the Alfven wave. They are 
shown by solid and dashed curves in panel a) respectively. The exact solution for v is shown 
by circles. In the slow and fast MHD waves variables p, u, B^, and p are all nonzero. The 
dimensionless pressure coincides in amplitude and phase with the density and is not shown 
in Figure [H Variables p, u, w, and i?^ are shown in panels b) and c) by solid, dash-dotted, 
dashed, and dotted curves respectively. The exact analytical solution for w is shown by 
circles. The analytical solutions for other variables coincide with the numerical curves and 
thus are not shown. We see that our numerical solutions reproduce the amplitudes, phases, 
and velocities of the Alfven, slow, and fast MHD waves very well. 



3. Results and discussion 

3.1. MHD waves generated by a single source in uniform and non-uniform 

background magnetic field 

In this section we present our results of numerical simulation of excitation and propaga- 
tion of MHD waves generated by a single source of vertical force with central frequency u = 
3.5 mHz placed at depth h^rc = 100 km in a rectangular region of size 15.5x15.5x12.5 Mm"^ 
(104x104x70 nodes). The horizontal grid is uniform with Ax = Ay = 150 km. The vertical 
grid is non-uniform. The grid step Az varies from 50 km near the photosphere to 600 km 
near the bottom of the computational domain. Time step At = 0.5 s was chosen to satisfy 
the Courant stability condition. Vector Bq = (Sq sin7, 0, cos7)^ of the uniform inclined 
background magnetic field lies in XZ plane and has inclination angle of 7 = 45° with respect 
to the top boundary normal. The magnetic field strength, Bq, varied from to 2500 G in 
our simulations. The ratio of the gas pressure to the magnetic pressure, plasma parameter 
f3 for typical sunspot penumbra values is shown in Fig. [2l 



- 8 - 



We describe here two examples of wave propagation in the uniform background magnetic 
field with a wave source located inside the magnetic region, and of a non-uniform magnetic 
field with the source located outside the magnetic regions. Since the convective motions 
on the Sun are inhibited in strong magnetic field regions, the acoustic sources there are 
suppressed. Thus, the second case better describes the realistic situation on the Sun. 

The simulation results for the uniform inclined magnetic field of Bq = 625 G are shown 
in Fig. [31 The vertical map of B'^ (panel c, right) reveals a strong Alfven wave, which 
propagates along the background inclined magnetic field lines. As expected, the Alfven wave 
is not presented in the map of density perturbations. The Alfven wave is generated due to 
the interaction between the wave source and the magnetic field at the source location. The 
concentric waves in the left panels represent a mixture of the fast MHD wave (analogous to 
p-modes in absence of the magnetic field) and the surface magnetic-gravity wave (analogous 
to f-mode). Since the wave speed depends on the angle between the vectors of magnetic field 
vector and wavenumber, the wave fronts are anisotropic. The separation of these two types 
of waves will be discussed in Sec. 3.2 together with analysis of phase and group travel times. 

The simulation results for the second case, when the source is placed outside the re- 
gion containing the background magnetic field are shown in Figure |H The domain size 
and grid parameters are the same, as in Figure [31 The background magnetic field Bq = 
(0, -Bo sin 7, Bq cos 7)"^ lies in YZ plane, has inclination angle 7 = 45° with the top boundary 
normal and is parallel to the planes showed in left columns of figure by dashed lines. The 
magnetic field Bq has constant values of and 2500 G in the regions on the left and right 
of the dashed lines respectively. Between the dashed lines the strength of the background 
magnetic field is linearly decreases from maximum value to zero. Such configuration of the 
background field satisfies condition (HyBq = 0. The magnetic field strength in this example 
is chosen higher than the typical penumbra value for a better demonstration of the magnetic 
effects. The wave source is placed in the region free from the background magnetic field. 
The waves generated by such a source are pure acoustic and surface gravity waves. When 
they enter the region occupied by the inclined magnetic field they are transformed into the 
fast MHD and slow magneto-gravity waves respectively. The Alfven wave does not appear 
in these simulations. Evidently, the fast MHD wave travels faster than the original acoustic 
wave, and its amplitude is reduced. In these simulations we do not notice significant wave 
transformation effects in the near surface reflection layers where the plasma parameter, /3, 
equal to 1. The helioseismic waves are trapped below the surface and according to our sim- 
ulations are not affected by transformation into other types (slow MHD and Alfven modes). 
We will discuss in more detail the role of the transformation for waves of different frequen- 
cies in a future publication. The main purpose of this paper is to discuss the effects of the 
inclined magnetic field on the observed travel time variations in the sunspot penumbra. 
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3.2. Phase and group travel time variations along the wavefront 

To study the travel time variations we performed simulations in a rectangular box of 
size 48x48x12.5 Mm^ (320x320x70 nodes) for different values Bq and inclination angles 
7 of the uniform background magnetic field: (625 G, 70°), (1400 G, 45°), and (1900 G, 
30°). The grid step size, source type and depth were chosen the same as in previous Section 
3.1. Total simulation time equals 6 hours of solar time. The simulation results for i?o=625 
G and 7=70° are shown in Figure O The top row represents k-u diagrams, the bottom 
one shows corresponded horizontal snapshots of the z-component of velocity at the level of 
the photosphere for the moment of time t = 30 min. The usual technique of fitting of the 
cross-covariance function by Gabor's wavelet was used to calcu late the travel times. This 



technique was developed for p-modes (IKosovichev fc Duvalllll997l ). The source of the vertical 



(z-component) of force generates a strong gravity wave (the lowest ridge of k-u diagram on 
panel a corresponds to the /-mode) which has to be filtered out. Results of separation of p- 
and /-modes are shown in panels (b) and (c) respectively. The maps of Vz for p- and /-modes 
(bottom row) were obtained by taking the inverse Fourier transform of the corresponding 3D 
spectra. It is clear, even without applying an /-mode filter, that starting form the distance 
of about 15 Mm from the source the p- and /-modes are spatially separated due to the 
different velocities. Due to the different dispersion relations the magnetic-gravity and fast 
MHD waves easily separated on time-distance diagram (see Figure E]). The solid black curve 
represents a theoretical time-distance curve for p-modes and standard solar model in absence 
of the magnetic field. We see, that the fast MHD waves are not significantly affected by the 
magnetic field a nd follow the theoretical p-m ode curve, calculated for the quiet Sun in the ray 



approximation (IKosovichev fc Duvalllll997l ) . The magnetic-gravity waves have characteristic 
"zebra" structure due to the difference between phase and group velocities. Their speed is 
less than the speed of the fast MHD wave. 

The magnetic field results in anisotropy of the wave properties. Therefore, the wave 
travel times measured from the line-of-sight component of the displacement velocity depend 
on the direction of the wave propagation and also on the viewing angle. We have investigate 
this effect for the case of the uniform inclined magnetic field. The choice of the coordinate 
system and geometry are shown in Figure [71 Horizontal XY-plane coincides with the pho- 
tosphere. The origin of the Cartesian coordinate system is placed at the point O above the 
wave source (the source itself is at the depth of 100 km below the photosphere). The uniform 
inclined background magnetic filed lies in the XZ-plane {By = 0) and has angle 7 = 45° with 
the normal to the photosphere. The location of point of observations Pobs is defined by the 
distance A from the wave source and the azimuthal angle a. The line-of-sight (LoS) direction 
is defined by two angles: angle 6 between LoS direction and local normal, and azimuthal 
angle ip between OX axis and projection of the LoS on the local horizontal plane. First, we 
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build 3D k-u diagrams for each velocity component u, v, and w, filter out the /-modes as 
shown in Figure [5] and calculate the Cartesian components of velocity perturbation Up, Vp, 
and Wp for p-modes by taking the inverse Fourier transform of corresponding spectra. Then, 
we calculate cross covariance 

C(n,r-2,r) = \ vtls{r,Mll{r,,t + T)dt (11) 

of LoS velocities 

v'flg = Up COS ipsinO + Vp sin tpsinO + Wp cos 6 (12) 

in the origin of the system of coordinates and in the observational point Pobs (A = 7.9 Mm). 
Then we fit the cross covariance function with Gabor's wavelet 

G(A, t) = A cos[ci;o(T — Tp)] exp 

where A = |ri — is the distance between points where LoS velocities v^Losi''^^^^) 
■^Losl^a,^) are measured, A is the amplitude, uq is the central frequency, Tp and Tg are the 
phase and group travel times respectively, and 6uj is the bandwidth. Parameters A, luq, Tp, 
Tg, and 6u! are free and have to be determined from the fitting procedure. Repeating this 
procedure for different observational points with the same A but different azimuthal angle 
we obtain Tp and Tg as functions of a. 

Travel time variations along the wave front for different Bq and different LoS angles 
are shown in Figure [Hi Panels a, b, and c correspond to Bq of 625 G, 1400 G, and 1900 G 
respectively. The solid, dashed, and dash-dotted curves correspond to = {0°, 90°, 180°} 
respectively. The angle between the LoS and the local normal is ^ = 20.5°. The mean phase 
travel times show variations of about 0.5 min along the wave front, the variation amplitude 
is a little smaller for strong (1400-1900 G) magnetic fields than for the weak 625 G field. 
For the weak magnetic field changes of ip change the shape of the curve, but not its average 
value. The strong magnetic fields cause anisotropy, and the mean value of the travel times 
changes with angle ip. Thus, averaging along the wave front for the strong magnetic fields 
gives variations of the observed mean travel times of about 0.1-0.3 min with angle ip (see 
detailed discussion in the next section). 

3.3. Azimuthal dependence of phase travel times in sunspots: comparison 

with observations 



(r - Tgf 



(13) 



Zhao fc KosovichevI (120061 ) observed different behavior of azimu thal dependence of ph ase 
travel times obtained from the SOHO/MDl Doppler shift data fIScherrer et al.l Il995l ) in 
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sunspots depending on their position on the disk. As far as we want to reproduce similar 
conditions (and geometry) in our simulations we give a detailed description of their algorithm 
of the phase travel time calculation. We will apply the same technique to our simulated data. 

The sunspot is located at latitude and longitude A on the eastern part of the disk as 
shown in Figure [H The sunspot meridian plane is shown by gray color. Orth i of the global 
system of coordinate with the origin in the center of the Sun is aimed at the center of the 
visible solar disk, orth k points to the north, and orth j points to the west. The local system 
of coordinates with the origin in the center of the sunspot is chosen in such a way that 
coincides with the local normal to the surface, is directed to the west, and Cy is directed 
along the meridian and forms a right vector triplet with and Cy. 

Observational point Pobs is chosen inside the sunspot penumbra. Azimuthal angle A is 
counted counterclockwise from local west direction e^- Two signals are calculated: the LoS 
velocity at Pobs and the LoS velocity averaged along the annulus (between inner and outer 
radii of a ring with average radius A = 8 Mm with the origin in the observational point). 
Cross covariance of these two signals is fit with the Gabors wavelet. Fitting procedure gives 
us mean phase and group travel times. We assume, that the annulus is small enough that 
sunspot magnetic field Bq inside the annulus can be considered as uniform. Hence the 
problem is reduced to simulations with the uniform inclined magnetic field shown in Figure 
m To compare results of the simulations with the observations we have to find from what 
direction we have to look at our simulation domain to match the observations. In other 
words, we have to find a relation between azimuthal angle A of the sunspot center and tp. 
Angle ip here has the same meaning as in FigureO This is the angle between the projection of 
the background magnetic field on the local horizon (XY-plane) and projection of the LoS on 
the same plane. The unit vector in LoS direction in the sunspot local system of coordinates 
is given by equation 

i = sin A — sin cos A + cos cos A e^. (14) 

The projection of * on local horizontal plane (defined by vectors e^. and Cy) has coordinates 
ip = (— sin A, — sin cos A, 0)^ (non-unit vector). Angle ip is given by equation 

i/j = —A — ZipE-j; = —A — arccos [ — j . (15) 

y V sin^ A + sin^ cos^ A J 

We performed numerical simulations of propagation of MHD waves in presence of the 
uniform inclined magnetic field for different inclination angles and strengths of the magnetic 
field. The goal of these simulations was to calculate contribution of the inclined magnetic 
field effect to variations of the mean travel times with the azimuthal angle. This is why we 
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used horizontally uniform standard solar model as a background model. The same model 
was used to obtain Figures [3l HI The mean travel times obtained for the photospheric level 
(panel a) and the level of 300 km above the photosphere (panel b) are shown in Figure [101 
The solid curve corresponds to the magnetic field Bq = 625 G with inclination angle 7 = 70°. 
The dashed and dash-dotted curves are corresponded to cases Bq = 1400 G, 7 = 45° and 
Bq = 1900 G, 7 = 30° respectively. The systematic shift of the calculated mean travel times 
with respect to the observations is caused by temperature effects (the sound speed is smaller 
inside sunspots than in the quiet Sun), which are not included in these simulations. The 
travel times systematically decrease when the strength of the magnetic field increases because 
the speed of the fast MHD wave increases with the magnetic field. For studying the inclined 
field effect we are mostly interested in relative variations of the mean travel times, provide 
the absolute values to give a general impression about the magnitude of contributions to the 
mean travel times from the temperature effect and variations of the fast MHD wave speed 
with the magnetic field. A model of a sunspot with both temperature and magnetic field 
effect will be presented in a future publication. 

The amplitude of the mean travel time variations for weak (625 G) magnetic field is 10 
times smaller than the observed quantity for the height of 300 km and even smaller for the 
photospheric level. For stronger magnetic fields (1400-1900 G) the amplitude of travel time 
variations is about 25% of the observed amplitude. Behavior of the travel times depends on 
the magnetic field strength and the level of observation of Doppler velocities. For the weak 
magnetic field (625 G) the phase of the travel time variations is opposite to the observations 
for both levels. For the height of 300 km above the photosphere where plasma parameter 
/3 <^ 1 simulated travel times show the same phase as in observations for both magnetic 
field strengths (1400 G and 1900 G). For the photospheric level and field strength of 1900 
G the phase of travel time variations coincides with the observations, while the travel time 
variations for Bq = 1400 G are in antiphase with the observations. Comparison of the curves 
at different heights for the same inclination angles and magnetic field strengths shows that 
for the photospheric layer the mean travel times are systematically smaller by only about 
0.12 min than for the layer of 300 km above the photosphere. Thus, the variations in the 
height of Doppler shift measurement do not have a significant effect. 

The curves in Figure [TD] were obtained for the annulus radius of 8 Mm and the sunspot 
located at heliospheric latitude = 19.5° and longitude A = —6.5°. The angle between LoS 
and local normal to the photosphere at the center of the sunspot {6 angle) is 20.5°. The 



parameters are chosen to match corresponded angles for Figure lb of IZhao fc Kosovichev 



fl2006[ ). 
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4. Conclusion 

The numerical 3D simulations of excitation and propagation of magneto-acoustic-gravity 
waves in the solar interior and atmosphere show that the presence of the inclined background 
magnetic field significantly alters properties of /-modes but has little effect on p-modes. The 
interaction of the wave source with the background magnetic field generates a strong Alfven 
wave. However, the Alfven wave does not appear when the wave source is located outside the 
magnetic region, and the acoustic-gravity waves propagate from the region without magnetic 
field into the region with magnetic field. 

The helioseismic travel times, obtained from cross covariance of the p-mode LoS veloc- 
ities at the observation point and the source point, show variations of about 1 min along 
the wave front (the amplitude depends on inclination 9 of the LoS). Due to the anisotropy, 
the travel time averaged along the wave front (like in the observational procedure) is not 
zero. Comparison of the variations of the mean travel times vs. the azimuthal angle of the 
observing point shows that the simulation results are in phase with the observations when 
Doppler velocities are taken at the level of 300 km above the photosphere (at the same height 
as the observed velocities are obtained). The travel time weakly depends on the height of 
observations. The amplitude of variations of the travel times obtained form simulations is 
about 25% of the observed amplitude even for strong fields of 1400-1990 G. It can be an 
indication that other effects (for example background flows or non-uniform distribution of 
magnetic field) can contribute to the observed travel time variations. The developed 3D 
MHD wave propagation code provides an important tool for further investigations of local 
helioseismology in regions with the strong magnetic field. 
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Fig. 1. — Comparison of analytic (cirves) and numerical (markers) solutions for plain Alfven 
(a), slow MHD (b), and fast MHD (c) waves respectively. For Alfven wave only v and By 
(marked by solid and dashed curves respectively) are non-zero. Circles mark an analytical 
solution for v. Solid, dash-dotted, dashed, and dotted curves in panels (b) and (c) correspond 
to p, n, and respectively. Analytical solution for w is shown by circles. 
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Fig. 2. — The ratio of the gas pressure to the magnetic pressure (plasma parameter (3) as a 
function of depth for the background model for three values of the magnetic field strength. 




(3) 






Fig. 3. — Snapshots of horizontal (left column) and vertical (right column) slices of the 3D 
domain for perturbations of density p (a), z- momentum PqW (b), and z-component of the 
magnetic field (c). Strong Alfven wave is generated due to the interaction of the wave 
source and the background magnetic field at the source location. In this example, Bo = 625 
G, and inclined by 45° in the XY-plane. 
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Fig. 4. — Snapshots of horizontal (left column) and vertical (right column) slices of the 3D 
domain for perturbations of density p (a), z- momentum PqW (b), and z-component of the 
magnetic field (c). The strength of the background magnetic field (inclined by 45° in 
the YZ-plane) is 2500 G in the region to the right from vertical dashed lines in left column 
and in the left region. The wave source is located in the region free of magnetic field. 
Propagating into the region with the background magnetic field acoustic and surface gravity 
waves are transformed into the fast MHD wave and magneto-gravity wave. The Alfven wave 
does not appear in this case. 
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Fig. 5. — Spectra {k-u diagrams) and corresponding maps of perturbation of the z- 
component of velocity are shown on top and bottom rows respectively. Results for the 
original MHD wave field and wave fields after filtering out /- and p-modes are shown on 
panels a, b, and c respectively. 
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Fig. 6. — Time- distance diagram for density perturbations. Solid curve represents a theoret- 
ical time-distance curve for p-modes and the standard solar model in absence of the magnetic 
field. Fast MHD wave and magnetic-gravity waves are separated due to the different disper- 
sion relations. 
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Fig. 7. — Choice of the coordinate system. The hne-of-sight direction is defined by angles ijj 
and 6. Position of observational point Pobs on the photosphere is fixed by azimuthal angle 
a and distance A from the projection of wave source O. 
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Fig. 8. — Mean phase travel times vs. azimuthal angle a calculated from Doppler LoS 
velocities at the hight of 300 km above the photosphere for different strengths of the uniform 
inclined background magnetic field (625 G for panel a, 1400 G for panel b, and 1900 G for 
panel c), and different LoS directions. Solid, dashed, and dash-dotted curves represent travel 
time variations for t/^ = 0°, 90°, and 180° respectively. 




Fig. 9. — To derivation of the relation between and A. Angle A is the azimuthal angle of 
the observation point in the local coordinate system associated with the sunspot. Angle ip 
is one of two angles determining the LoS direction. 
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Fig. 10. — Phase travel time obtained from simulations averaged along with annulus with 
radius 8 Mm as a function of azimuthal angle A. The solid, dashed, and dot-dashed curves 
corresponds to the following combinations of the magnetic field strength and inclination 
angles: (i) Bq = 625 G, 7 = 70°, (ii) Bq = 1400 G, 7 = 45°, (iii) Bq = 1900 G, 7 = 30°. 
Data points with errorbars represent observations. 



